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PARALLEL  AND  INTERSECTING  FLATS  FRACTIONS: 


EST I MAB I L I TV  AND  ALIAS  STRUCTURES 

Donald  A.  Anderson 
and 

Ann  M.  Thomas 


ABSTRACT 

Regular  fractional  factorial  designs  for  the  sn,  s = pa. 


factor  ial-^may  be  obtained  as  solutions  sets  to  linear  equations, 
of  the  form  A £ = £ over  GF(s).  ’“'It  is  well  known  how  to  deter- 
mine which  factorial  effects  are  estimable,  and  how  to  construct 

Sys  Lc  iy\ 

the  alias  sets  directly  from  the^matr ix;A.  Consider  the  k sets 
of  equations  over  GF(s). 


A • t_  = c . , i = l ,2 k 

k 

where  A.  is  m.  x n of  rank  m. , and  let  T = .U,{t  A.t  = c.}. 

If  Aj  = A^  = ...=  A^,  T is  called  a parallel  flats  fraction,  and 
if  not  all  A.  are  equal,  T is  called  an  intersecting  flats 
fraction.  ^There  is  no  unified  theory  developed  for  determining 
estimability  and  alias  structures  from  the-se- types  of  fractions,.-  irvoW'r 

q Fkts.  r 

'■The  purpose  of  this  paper  is  to  present  some  preliminary 
results  in  the  development  of  such  a theory. 


1 . Fractional  Replication 

As  the  number  of  factors  in  an  s'1  experiment  increases,  the  high 
cost  of  experimentation,  expecially  in  industry,  often  prohibits  use  of 
a complete  factorial  design.  If  economy  of  time,  space,  and  materials 
is  an  important  consideration,  and  if  it  is  reasonable  to  assume  that 
certain  higher  order  interactions  are  negligible,  then  a fractional 
factorial  design  may  be  appropriate.  The  concept  of  fractional  repli- 
cation for  2n  and  3n  experiments  was  introduced  by  Finney  (1945)  as  a 
logical  outgrowth  of  the  theory  of  confounding,  and  since  that  time  the 
theory  of  fractional  replication  has  developed  rapidly.  Experimenters 
involved  in  high  cost  research  have  found  the  designs  to  be  useful, 
particularly  as  screening  experiments  run  on  a large  number  of  factors. 

The  notion  of  linear  flats  within  a finite  Euclidean  geometry,  as 
developed  by  Bose  and  Kishen  (1940)  and  Bose  (1947),  is  a useful  device 
for  partitioning  degrees  of  freedom  and  constructing  fractional  factorial 
designs.  Suppose  s = p , where  p is  prime,  and  consider  an  s factorial 
experiment  involving  n factors,  each  at  s levels.  Let  every  possible 
treatment  combination  or  "run"  be  represented  as  an  n x 1 vector  t 
with  elements  from  a Galois  field  of  order  s,  denoted  GF(s),  and  with 
ith  element  corresponding  to  the  level  of  the  ith  factor.  Then  there 
clearly  exists  a one-to-one  correspondence  between  the  s"  possible 
treatment  combinations  and  the  sn  points  in  the  Euclidean  geometry  of 
order  n over  the  field  of  order  s,  denoted  as 

EG(n,  s)  = {(tj,  ty,  ...,  tn)  | tj  €.  GF(  s ) , i = 1,  ...,  n). 


2 

A linear  (n  - k)-flat  in  FG(n,  s)  is  defined  to  be  the  set  of  points  in 

FG(n,  s)  that  satisfy  a consistent  system  of  k independent  linear 

n 

equations,  each  of  the  form  J a.t.  = a , where  a-  = 0,  a, , a 

1 l i r’  0 ’ 1’  ’ s-1 

i=l 

denote  the  elements  of  GF(s)  and  a = a.  for  some  i = 0,  1,  ors-L 

r i ’ ’ 

Consider  a pencil  of  s parallel  (n  - l)-flats  with  equations  of  the 
n 

form  } a.t.  = a , where  a.  is  fixed,  i = 1,  ...,  n,  and  a assumes 
. L.  l l r’  l ’ ’ ’ ’ r 

1 = 1 

consecutively  the  s values  a(),  a^ , ...,  a j.  Any  such  pencil  partitions 
the  sn  points  in  EG(n,  s)  into  s subsets  of  sn  ^ points  each,  with  con- 
trasts between  these  subsets  accounting  for  s-1  degrees  of  freedom. 
Thus  the  pencil  t^  = a , i fixed,  and  r = 0,  1,  . ..,  s - 1,  is  used  to 
define  the  (s-1)  degrees  of  freedom  corresponding  to  the  main  effects 
of  the  ith  factor,  the  pencils  given  by  equations  of  the  form 

tj  + aut  i ' = ar  » 1 an^  i ' fixed , u = 1 , 2,  ...,s-l  and  r = 0,  1,  ...,  s - 1, 

2 

are  used  to  define  the  (s-1)  degrees  of  freedom  corresponding  to  the 

two-way  interaction  of  Factors  i and  i',  and  so  on.  Suppose,  for 

instance,  that  each  factor  in  a 3n  experiment  appears  at  levels  0,  1 , and 

2.  The  two  degrees  of  freedom  provided  for  the  main  effects  of  the  ith 

factor  by  the  pencil  t.  = r,  r = 0,  1,  2,  correspond  to  the  linear 

contrast,  written  as  (t  | t^  5 2}  - (t  | t,  = 0),  and  to  the  quadratic 

contrast,  written  as  {t  I t.  = 2}  - 2 { t I t.  = 1}  + {t  I t.  = 0}.  For 

any  two  factors  F^  and  F^,  the  same  type  of  contrasts  may  lie  used  to 

partition  the  two  degrees  of  freedom  for  interaction  F^F^,  specified  by 

2 

the  pencil  t^+t^,  =r,  r=0,  1,  2,  and  the  two  degrees  for  F^F^ , 
specified  by  the  pencil  t..  + 2t^,  = r,  r = 0,  1,  2.  This  procedure  gives 
a single  degree  of  freedom,  breakdown  of  the  four  degrees  of  freedom 
corresponding  to  the  interaction  of  Factors  i and  i'.  Clearly,  the 


f 


3 


components  of  this  breakdown  are  not  the  same  as  the  components  of  the 
breakdown  into  single  degrees  of  freedom  for  linear  x linear,  linear  x 
quadratic,  and  so  on,  that  was  described  in  Section  1.1. 

The  simplest  kind  of  fractional  factorial  design  for  the  sn  experi- 
ment consists  of  all  treatment  combinations  corresponding  to  points  in 
EG(n,  s)  that  lie  on  a particular  (n  - l)-flat  in  a pencil  of  s parallel 
flats.  The  aliasing  structure  for  such  a design  is  easily  determined  by 

taking  linear  combinations  of  the  defining  relationship  for  the  selected 
n 

pencil,  say  £ a.t.  = a , r = 0,  1 , . . . , s - 1 , with  the  defining 
i=l  1 1 r 

relationships  for  main  effects  (t^  = ap,  r = 0,  1,  ...,  s - 1),  two- 
factor  interactions  (t.+at.,=a,r=0,l,  . ..,  s - 1),  and  so  on. 
In  general,  the  defining  pencil  should  be  selected  so  that  effects  of 
interest  are  aliased,  or  confounded,  with  negligible  effects. 

4 

For  example,  one  fraction  of  a 3 experiment  with  Factors  F^ , F^ , F^, 
and  F^,  each  at  levels  0,  1,  and  2,  consists  of  the  27  points  on  the 
three-flat  specified  by  tj  + 2t?  + t3  + 2t^  = 0.  Equivalently,  this 
fraction  is  often  represented  by  the  defining,  relations 


I 


o o 99 

F F F F (=  F F F F ). 
1 2 3 4 1 2 3 4 


2 

In  this  design  the  F F„  interaction,  winch  by  definition  is  based  on 

£ #0 

contrasts  among  sums  of  points  for  which  t^  + 2t^  = constant,  is  aliased 
2 

with  F^F^  since  (t^  + 2t^  + t3  + 2t4)  + (t^  + 2tg)  = t^  + 21^  = constant. 
Similarly,  the  F interaction  is  aliased  with  the  F^F^  interaction 

since  (t^  + 2t^  + tg  + 2t4>  + 2(1^  + tj)  = 2(t^  + ^ three-factor 


and  higher  order  interactions  are  negligible,  the  alias  sets  for  this 


4 


fractional  (m),  {F^,  (F^,  (F^,  {F^ } , {F1F2},  (Firi4  > , {F^},  {F^}, 

('VV  r2r„),  (FjFg},  fr2r*).  tqr’,  ly*),  (r,^,  r2r*l. 

In  general,  regular  sn  ^ fractions  of  the  s11  factorial  experiment 
consist  of  solutions  over  GF(s)  to 


At  = c 


(1) 


where  A is  a k x n matrix  of  rank  k and  c_  is  k x 1,  both  over  a Galois 

field  of  order  s.  Geometrically,  the  equations  (1)  represent  a linear 

n k 71  n 

(n  - k)-flat  of  s points  in  FG(n,  s).  If  7 a.t.  = a and  7 alt. 

,L , l l r . la  i i 

i=l  i=l 

= a 1 are  two  equations  used  to  define  a fraction,  their  generalized 

n n 

interaction  is  specified  by  ( ) a.t.)  + a ( 7 a',  t.)  = a , u = 1,  ...,  s - 1 , 

1 i i u ,L . 11  w 

i=l  i=l 

w = 0,  ...,  s - 1,  and  it  is  an  easy  matter  to  determine  the  alias  sets 

for  the  fraction  by  considering  all  linear  combinations  involving 

defining  relationships  and  their  generalized  interactions. 

n ” k 

In  practice  a regular  s fraction  is  often  specified  by  defining 

relations  in  k "words",  chosen  so  that  effects  of  interest  are  aliased 

G ”■  2 

with  higher  order,  negligible  interaction  effects.  Fon  example,  a 3 

4 

fraction  consisting  of  the  N = 3 =81  treatment  combinations  t that  are 

solutions  to 


f1  1 

1 

1 

0 

°lt  =frll 

o 

o 

1 

1 

2 

1 

ij-  lc2l 

can  be  specified  by  the  defining  relations 


? 2 

r r p p = F F F F = F F F F r 
1234  3 4*5  6 12356’ 


L 


V 
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where  * 6 ^le  Senera^^zec^  interaction  generated  by  F ^ i 2r  3^  4 

and  F F!T  F . The  resulting  fraction  is  a design  in  which  main  effects 
3 4 D o 

are  aliased  with  three- fact or  and  higher  order  interactions.  Thus,  il 
three-factor  and  higher  order  interactions  are  negligible,  the  design 
provides  estimates  of  main  effects  in  the  presence  of  non-negligible 
two-factor  interactions. 


Bose  (1944)  defined  a linear  combination  3 of  the  parameters  in  a 
linear  model  to  be  estimable  if  and  only  if  there  exists  a linear  com- 
bination c' Y of  the  observations  such  that  F(£'X^  = A'J1‘  A complete 
factorial  design  provides  an  estimate  for  every  single  degree  of  freedom 
component  of  6.  With  a fractional  factorial  design,  on  the  other  hand, 

B.  s B_  is  estimable  if  and  only  if  the  column  of  X that  is  multiplied  by 
B^  in  tea  model  E(Y_)  = X_3  is  linearly  independent  of  all  remaining 
columns  of  X. 

n-}< 

Box  and  Hunter  (1961),  in  a paper  dealing  with  regular  ? fractions 
of  the  2n  factorial,  introduced  the  term  resolution  of  a design  as  a 
means  of  classifying  fractional  factorial  designs  on  the  basis  of 

n ■”  k 

estimability  of  effects.  The  term  applies  directly  to  the  regular  s 
fractions  of  the  sn  factorial  experiment,  in  which  case  the  design  is 
resolution  r if  all  linear  combinations  of  the  rows  of  A from  equation 
(1)  have  at  least  r nonzero  coordinates.  The  names  resolution  III, 
resolution  IV,  resolution  V,  and  so  on,  make  obvious  reference  to  the 
particular  kind  of  fractions  considered  by  Box  and  Hunter.  For  example, 
a resolution  IV  design  for  the  3°  factorial  experiment  requires  each 


6 


word  in  the  defining  relations  to  have  at  least  four  letters  and  at  least 
one  word  to  have  exactly  four  letters.  Such  a design  permits  estimation 
of  the  general  mean  and  all  main  effects  if  all  three-factor  and  higher 
order  interactions  are  negligible.  Further,  the  estimates  are  uncorre- 
lated with  each  other. 

To  the  extent  that  the  name  resolution  r still  carries  implication 
to  the  regular  sn  ^ fractions  of  the  s"  experiment  it  is  an  unfortunate 
terminology.  A slightly  modified  definition  which  is  now  generally 
accepted  was  given  by  Webb  (1965,  1968)  and  Margolin  (1969a)  as  follows. 
Definition  1.  If  a design  is  such  that  all  effects  (main  effects, 
interactions)  involving  r or  fewer  factors  are  estimable,  ignoring  all 
interactions  of  r + 1 or  more  factors,  the  design  is  said  to  be  of 
resolution  2r  + 1 •,  if  all  effects  involving  r - 1 or  fewer  factors  are 
estimable,  ignoring  all  interactions  of  r + 1 or  more  factors,  the  design 
is  said  to  be  of  resolution  2r. 

Note  that  this  definition  refers  only  to  the  estimability  of  effects , 
not  to  any  specific  method  of  construction.  Further,  it  does  not  require 
estimates  to  be  uncorrelated  with  each  other  or  u to  be  estimable.  In 
general,  designs  of  odd  resolution  permit  estimation  of  all  effects  not 
assumed  to  be  zero,  wh,ile  designs  of  even  resolution  permit  estimation 
of  certain  effects  in  the  presence  of  other  non-zero,  nonest imable 
effects.  In  practice,  designs  of  resolutions  III,  IV,  and  V are  of  the 
most  interest.  A resolution  III  design  allows  estimation  of  main  effects 
when  two-factor  and  higher  order  interactions  are  negligible,  and  a 
resolution  V design  allows  estimation  of  main  effects  and  two-factor 
interactions  when  three-factor  and  higher  order  interactions  are 


negligible.  A resolution  IV  design,  on  the  other  hand,  permits  esti- 
mation of  main  effects  in  the  presence  of  nonestimable  two-factor 
interactions  when  three-factor  and  higher  order  interactions  are 
negligible. 

This  paper  was  motivated  by  the  authors'  search  for  resolution  IV 
designs  for  t lie  Sn  factorial  (Anderson  and  Thomas  1975a,  b).  This  search 
led  to  the  consideration  of  parallel  and  intersecting  flats  construction. 
At  the  present  time  there  is  no  unified  theory  developed  for  determining 
the  alias  sets,  or  even  which  effects  are  estimable,  from  such  con- 
structions. The  purpose  of  this  paper  is  to  begin  a development  of 
estimability  and  alias  structure  theory  for  parallel  and  intersecting 
flat  fractions.  Because  of  the  initial  motivation  most  of  the  examples 
relate  to  estimability  of  main  effects  in  the  presence  of  two-factor 
interactions . 

2 . Construction  Using  Intersecting  Flats 

In  Section  1 a regular  s'1  ^ fraction  of  the  sn  factorial  experiment 
was  defined  to  be  the  set  of  solutions  over  GF(s)  to 

At  = c,  (3) 


where  A is  a k x n matrix  of  rank  k and  c is  k x 1,  both  over  a field  of 


n-k 

order  s.  Thus  there  is  a one-to-one  correspondence  between  the  s 
points  of  EG(n,s)  that  lie  on  the  linear  (n  - k)-flat  specified  by  (2) 

n-k 

and  the  s treatment  combinations  that  are  in  the  fraction  determined 
by  (2).  An  alternative  to  choosing  points  on  a single  flat  in  EG(n,s) 
is  to  take  the  union  of  thepoints  on  several  flats.  Thus  consider  the 
flats  generated  by  equations 


A.t  = £i>  i = 1,2, . . . ,r  (3) 

where  A.  is  m.  x n of  rank  m.  and  c.  is  m.  x 1.  The  design  T corresponding 
li  1-11 

to  (3)  is 

r 

T = U { 1 1 A. t = c.}  . (4) 

i=l  - -1 

The  i^-h  flat  contains  sn  m-*-  points,  but  since  the  flats  may  intersect  in 
various  ways  the  number  of  points  in  T as  well  as  the  estimability  of 
factorial  effects  depend  on  the  A^  and  c^  in  a rather  complex  manner. 

Consideration  of  designs  of  type  (4)  was  motivated  by  a search 
for  a general  series  of  minimal  or  near  minimal  resolution  IV  designs  for 
the  sn  factorial.  Theorems  1 and  2 produce  such  a series  in  N = s(s  - l)n 
runs,  only  s(s  - 2)  more  than  the  perhaps  unattainable  lower  bound. 


Theorem  1.  Let  A,,  A„ A be  (n  - 2)  x n matrices  of  rank  (n-2) 

12  n 

such  that  for  1=1,  2,  ...,  n, 

1.  the  iTh  column  of  A.  is  0, 

l — 

2.  one  column  of  A^  has  all  nonzero  elements,  and 

3.  the  remaining  (n  - 2)  columns  of  are  some  permutation  of  I^n_2) 


9 


Then  the  fraction 


n 

T = U ti  I A;i  = £• 1 (5) 

i = l 

is  resolution  IV  for  any  choice  of  the  c^. 

froof.  Each  row  of  A.  contains  exactly  two  nonzero  elements  and  has  ith 

l 

element  zero.  Further,  any  linear  combination  of  two  or  more  rows  of  A. 
will  have  at  least  two  nonzero  elements  and  a zero  in  the  ith  position. 
It  follows  directly  that  in  the  ith  set  of  equations,  the  main  effect  of 
the  ith  factor  is  aliased  only  with  three  factor  or  higher  order  inter- 
actions. Thus  the  main  effects  of  the  ith  factor  are  estimable  from  the 
runs  corresponding  to 


A.t_  = £.  , i = 1,  2,  ...,  n. 


so  all  main  effects  are  estimable  from  the  union 


n 

T = U U 
i=l 


A . t 

l— 


c . } . 

—l 


The  proof  is  complete. 


Example  j . The  design  T formed  by  taking  the  union  of  solutions  over 

5 

GF(s)  to  the  following  five  systems  of  equations  is  an  s resolution  IV 
design  in  s(s  - 1)5  runs. 


Flat  1 . Flat  2.  Flat  t. 


9 

-1 

1 

n 

o' 

o’ 

’-1 

0 

1 

0 

*1 

0 

V 

-1 

1 

g 

0 

1 

0 

'o 

0 

-1 

0 

l 

0 

t = 

0 

-1 

0 

0 

1 

0 

t = 

1 

-1 

0 

0 

.1 

0 

t:  = 

1 

0 

-1 

0 

0 

Pi 

pj 

-1 

0 

0 

0 

1 

1 

-1 

0 

0 

0 

p 

1 

Flat  4.  Flat  5. 


-1  1 0 0 o' 

"o' 

-1  1 0 0 o' 

0 

-10100 

1 = 

0 

-1010' 

t = 

0 

-10001 

1 

-1  0010 

.0.. 

iO 


The  union  of  the  points  on  these  five  flats  with  s = 3 is 

012  01?  012  012  012  012  012  012  012  01?' 

01?  201  120  201  012  01?  012  012  012  01 2 

T = 012  201  120  120  120  201  012  012  01?  012 

012  201  120  120  120  120  120  201  012  01? 

012  201  120  120  120  120  120  120  120  201 


where  the  underscoring  indicates  the  treatment  combinations  generated  by 
the  ith  system  of  equations,  i = 1,  2,  3,  4,  5. 

Theorem  2.  For  s = p , p prime,  there  exist  fractions  of  resolution  IV 
in  N = s(s  - l)n  runs  for  the  s11  factorial. 

Proof . The  proof  is  by  construction.  Let  the  column  of  all  nonzero 

elements  of  the  in  Theorem  1 be  -Jj_,  that  is,  an  (n  - 2)  column  with 

every  element  -1  e GF(s).  Further,  let  the  column  -1  be  tiie  second 

column  in  and  the  first  in  A0,  A^,  ...,  An-  Finally,  arrange  the 

rows  of  each  A.  so  that  the  set  of  (n  - 2)  columns,  each  with  a single 

one  and  (n  - 3)  zeros,  in  order  form  I . The  arrangement  in  Example 

n-  2 

3.1  illustrates  this  particular  ordering.  The  c_. , i = 1,  ...,  n,  are 
then  selected  so  that  flats  1 and  n intersect  in  s points,  flats  i and 
(i  + 1),  i - ...,  n - 1,  intersect  in  s points,  and  all  remaining 

intersections  of  pairs  of  flats  are  empty.  The  construction  is  repre- 
sented pictorially  in  Figure  1. 


i. 
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ii 


n.u  1 


Flat  (i  + 1) 

Figure  1.  Geometric  Representation  of  Construction 
for  Design  in  s(s  - l)n  Runs 

With  c_j.  = (c^.,  c2i’  cn_2  i^'’  * = •••’  n ’ t^le  construction 

is  achieved  as  follows. 

1 * £2  “ for  some  constant  a implies  flats  1 and  2 inter- 

sect in  s points. 

2'  Ck+l,n  = cki  + c j n ’ k = 1 ’ . ..,  n - 3 implies  flats  1 and  n 
intersect  in  s points. 

3-  cki  = ck>i+i‘>  k = l,  ...,i-2,  i,  ...,n-2,  and  i = 3, 

n - 1,  implies  flats  i and  (i  + 1)  intersect  in  s points,  i = 1, 

. . . , n - 1 . 

4-  cn_2,i  ^ cn-2  |+cjj’a=3’---»n-1  implies  flat  1 does  not 
intersect  with  flat  i , i = 3,  . . . , n - 1 . 


5.  c„  , = c.  ^ . , 1 = 2,  . . . , n - 2,  and  i'  = i + b,  b > 1, 

implies  flat  i does  not  intersect  with  flat  i + b,  b >1,  i 

n - 2. 

For  arbitrary  £1  = [c^  , c ,]  and  , and  any  c.^  t c1?  - c 

the  following  construction  clearly  satisfies  conditions  1-3. 

ck2  C12  + ckl  ' cll’  k " 3’  •••■>  n - 2 

cli  = Cln’  1 = 3’  •••’  n " 1 

Cki  = Ck,i-1’  k = 1 - n ' 2’  Cln  + Ck-l,l’  K = 2>  1 

i.  = 3 , . . . , n . 


Also,  since  cln  i c 2 - c ^ implies  c1n  i c^_0  0 - c 


In  n-2,2  n-2,1 


implies  c. . i c 


1 i n-2,1  n-2, 1 


, i = 3 , . . . , n- 


implies  c „ . / c _ , + c.  . , i = 3,  ... 
n-2,i  n-2,1  li 

n - 1, 


condition  4 is  satisfied.  And  for  planes  i and  i'  = i + b,  b > 1,  i : 

n-2’  °ln^c12-cll  imPlies  ci-l,l  + Cln  * Ci-l,l  + (C12  * 
implies  cln  + c._1?1  i c._l5? 


implies  c . . . i c . , . 

ii'  i-l,i. 


so  condition  5 is  satisfied. 

One  such  choice  of  the  c.  is  c.  = c = 0,  c_  = 1,  c0  = (0,  1,  1, 

— i — i — n — — Z — — o 

...,  1)’,  c.  = (0,  0,  1,  ...,  1)’,  ...,  c = (0,  0,  0 1)'. 

— — n - 1 

Example  2.  The  3^  design  in  Example  1 illustrates  the  construction  of 
Theorem  2. 

The  designs  constructed  by  Theorem  2 have  a convenient  represents 
tion  in  terms  of  parallel  one-flats.  Since  flats  1 and  n,  as  well  as 
flats  i and  i+1,  i=l,2, . . . ,n-l,  intersect  in  s points,  each  of  the 


13 


! 


system  of  equations 


A, 

fc,! 

r a . i 

c . 

1 

t = 

-l 

, and 

i 

t = 

-1 

A 

n 

C 

-n_ 

*i+l 

i = 1,  2,  . . . , n-1,  (6) 


is  of  rank  n-1  and  consistent.  In  fact,  each  can  be  reduced  to  the 
form 

t = d . (7) 

Ihe  equation  (7)  with  the  (i  - l)th  row  eliminated,  i = 2,  3,  n, 

are  precisely 

Ait  = d*  , 

where  d*  is  the  (n  - 2)  * 1 vector  obtained  from  d by  eliminating  the 
(i  - l)th  element.  In  relation  to  Theorem  2 d*  = c..  Thus  the  set  of 

- -l 

points  on  the  ith,  i = 2,  ....  n,  flat  may  be  represented  as  s parallel 
one-flats  defined  by 

t = d^,  k = 0,  1,  2 s-1,  (8) 

where  the  (i  - l)th  element  of  d^  is  e GF(s)  and  the  remaining  n - 2 
elements  are  d*  = c^.  Similarly,  if  the  first  equation  of  (7)  is  subtracted 
from  each  of  the  remaining,  the  result  is 


14 


Al^  = 


d2  -dl 
d3  “dl 


d - d 
n-1  1 


= si  • 


(9) 


and  the  s points  on  A^_t  = are  obtained  by  keeping  fixed  and  let- 
ting d^  take  all  values  in  GF(s).  Thus  in  the  parallel  flats  representa- 
tion (7)  the  coefficient  matrix  is  fixed  while  vector  d is  varied. 


Example  3.  The  parallel  flats  representation  of  the  y design  of  Example 
1 is  given  by 


-11000 

-10100 

-10010 

-10001 


t = d^ , i = 1,  2 10 


where  the  d.  are  the  columns  of 

-l 


D = 


0 

2 

1 

2 

0 

0 

0 

0 

0 

0 

0 

2 

1 

1 

1 

2 

0 

0 

0 

0 

0 

2 

1 

1 

1 

1 

1 

2 

0 

0 

0 

2 

1 

1 

1 

1 

1 

1 

1 

2 

The  one-flats  of  each  plane  are  identified  by  the  spacing  in  the  array 
T of  Example  1.  Several  additional  examples  are  included  in  Section  3. 

4 

3.  The  3 Experiment  - More  Examples. 

4 

A somewhat  detailed  consideration  of  the  3 experiment  will 
illustrate  both  the  versatility  and  the  complexity  of  design  construct  ion 
from  intersecting  linear  flats. 


t- 


— 
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For  the  3 factorial  consider  first 


A1  " 


0 2 10 
0 2 0 1 


2 0 10 
2 0 0 1 


» Ay 


2 10  0 
2 0 0 1 


2 10  0 
2 0 10 


(10) 


as  specified  in  Theorems  2 and  3 of  Section  2, 


-i  = (cli’  c2i)'>  1 = 1 > 2 » 3, 


and 


T = U { 1 1 A . t = c , } . 
i=l  " i_  _1 


The  number  of  runs  in  T,  which  corresponds  to  the  number  of  points  in  the 
geometric  representations  of  T,  clearly  depends  on  the  choice  of  c.,  i = 

1,  2,  3,  4.  For  instance,  if  c.  = 0,  i = 1,  2,  3,  4,  the  design  includes 
the  27  runs  represented  in  Figure  2 of  Section  3.  On  the  other  hand,  if 
the  c^  are  chosen  to  be 
c’  = (0,0) 
c’  = (1,1) 
c’  = (0, 1) 

£4  = (0,0) 

the  result  is  a design  in  24  runs.  In  general,  with  , Ayt  Ay,  and 

defined  as  above, any  choice  of  yields  a two-flat  which  is  simply  the 

union  of  three  non-intersecting  one-flats,  each  being  of  the  form 

4 

{t,  t + 1,  t + (2)1}.  Thus  a design  T = U { t | A . t = c.},  expressed  here 

i=i  ~ 1- 

as  the  union  of  four  intersecting  or  nonintersecting  two-flats,  could 
alternatively  be  written  as  a union  of  parallel  one-flats.  That  is, 

4 r 

T = U { 1 1 A t = c } = U { 1 1 Bt  = d } , 

i=l  " i=l  " 


10 


where  B is  a 3 x 4 matrix  of  rank  3 and  r depend';  on  the  inter:, eot ion 
structure  of  the  two-flats.  As  in  Section  2 , it  is  convenient  to  let 


B 


-110  0 
-10  10 
-1001  . 

_ 


(11) 


The  examples  which  follow  further  illustrate  t lie  dependence  of  the 

number  of  points  in  T on  the  choice  of  c.  , e,,,  c„,  and  c,  . for  each 

— I — / — j — 4 

example  it  is  assumed  that  A.,  i = 1,  2,  3,  4,  and  B are  as  defined  in 
( 10  and  ( 11  ) , respectively.  Thus  C = [c^  c . ^ ] completely 

specifies  a design  as  a union  of  two-flats,  while  D = [d^  d_j  •••  d.r3 
specifies  it  as  a union  of  one-flats.  In  the  geometric  configurations 
lines  represent  two-flats  and  points  represent  one-flats. 


L'xample  3. 


, . _ fo  0 1 o' 

36  run  design  c = 0 2 l 


D 


012012111000 
- 012000012111 
201000222012 


Figure  2 ■ 


Geometric  Representation  of  36  Run 


Dos i gn 


Example  4 • 


33  run  design  0 


D 


01212111000 
= 01200012111 
01200222012 


Figure  3.  Geometric  Representation  of  33  Run  3 Design 


Figure  5.  Geometric  Representation  of  27  Run  3 Design 


Ifc 


Clearly,  this  27  run  design  is  not  the  same  as  the  27  run  design 
pictured  in  Figure  2 of  Section  3.1. 

There  is  no  uniqueness  to  the  particular  constructions  presented  in 
Theorems  one  and  two  of  Section  2.  Other  choices  of  1 he  A.  inav 

i 

also  give  rise  to  resolution  IV  designs.  For  example,  the  27  run  3* 
foldover  designs,  [given  by  Anderson  and  Thomas  (!9/5a)],  all 
have  representations  as  four  intersecting  two-flats.  The  A.  , i = l,  2,  3,  4, 
for  these  designs  are  listed  in  Example  7.  In  each  case  - c3 

= S,  = 0. 

4 

Example  7.  A.  matrices  for  27  run  3 resolution  IV  Foldover  Designs 


(a)  Design  2 (orthogonal  array  of  strength  3) 


(b)  Design  3 


(c)  Design  4 


The  smallest  known  designs  for  t he  3 experiment  are  the  ones  gener- 
ated by  Theorem  two  of  Section  2 • Since  these  designs  require  N - 24 

. 4 

runs  while  Margolin's  lower  bound  for  resolution  IV  3 designs  is  N > 21, 
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it  is  of  interest  to  investigate  the  possibility  of  using  four  intersec- 
ting two-flats  to  construct  a 21  run  resolution  IV  design. 

Consider  the  21 -point  configuration  pictured  in  Figure  6. 


Figure  6-  21-point  Configuration  in  CG(4,  3). 


In  Figure  6 every  pair  of  two-flats  intersect  in  exactly  three  points, 
and  every  triple  of  two-flats,  except  the  one  consisting  of  flats  1,  2, 
and  3,  intersect  in  exactly  one  point.  The  intersection  of  flats  1,  2, 
and  3 is  empty.  For  i = 1,  2,  3,  and  4 the  ith  flat  consists  of 
solutions  over  GF(3)  to 

Ajt_  = c.  (12) 

where  A.  is  2 x 4 of  rank  two  and  c.  is  2 x 1 . Thus  any  two  of  the  four 

l —l 

flats,  say  flats  i and  j , intersect  in  exactly  t liree  points  if 


rank 


rank 


3, 


(13) 


and  any  three  of  the  four  flats,  say  flats  i.  i,  and  k,  intersect  in 


20 


exactly  one  point  if 


A. 

A. 

c . 

l 

l 

—1 

A . 

= rank 

A. 

c . 

) 

I 

-J 

Ak 

\ 

Flats  1,  2,  and  3 have  an  empty  intersection  if  and  only  if 


IV 

r 

Ai 

s-1 

rank 

A2 

= 3 and  rank 

A2 

-2 

_A3_ 

_A3 

£3 

(14) 


(15) 


In  accordance  with  Theorem  1 of  Section  2 let  A.,  i = 1,  2,  3,  4, 
be  2 x 4 matrices  of  rank  2 such  that  the  ith  column  of  A.  is  0,  one 

l — 

column  of  A^,  denoted  by  cu  = (a^.,  a^,-)',  has  all  nonzero  elements,  and 

the  remaining  two  columns  of  A.  form  I-.  Since  Theorem  1 guarantees 

4 

that  a fraction  T = U {t_  | A . t_  = c_.  } constructed  using  these  A.  as 

i = l 1 1 1 

coefficient  matrices  is  resolution  IV,  it  is  of  primary  interest  to  know 

whether  for  i = 1,  2,  3,  4 there  is  some  choice  of  a.  = ^ali’  a2i^'  dnt^ 

— i = ^cli’  C2i^'  ^rom  GF(3)  that  will  yield  the  construction  of  Figure  6. 

The  result,  of  course,  would  be  a minimal  resolution  IV  design  for  the 

4 

3 experiment.  Unfortunately,  if  the  A^s  are  defined  as  in  Theorem  1 
if  condition  (13)  is  satisfied  for  every  pair  of  two-flats,  and  if 
condition  (14)  is  satisfied  for  flats  i,  j , and  k when  (i,  j,  k)  = 

(1,  2,  4),  (1,  3,  4),  and  (2,  3,  4),  then  it  is  impossible  for  flats 
1,  2,  and  3 to  satisfy  condition  (15). 

Since  the  A.s  specified  by  Theorem  1 cannot  be  used  to  produce 
the  configuration  of  Figure  6,  consider  instead 


r: 

1 0 

a.  • 

b,  .1 

= 

ll 

li 

Li 

3 1 

a2i  b2iJ 

(lb) 


21 


and  the  corresponding  two-flats  in  CG(4,  3)  defined  by  A^t_  = c^ , 
i = 1,  2,  3,  4.  With  the  A^s  defined  in  this  way  it  is  possible  to 
select  from  Gf(3)  values  for  a.  = (a^,  a^.)',  b.  = (b,,  . , bo;  )',  and 


2 1 ' 


— i = ^Cli’  C2i^'’  i = 1>  2,  3,  4,  so  that  the  resulting  two-flats 
intersect  as  in  figure  6.  However,  an  exhaustive  consideration  of  all 
possibilities  yielded  only  two  intrinsically  different  constructions  of 
this  kind,  the  conditions  for  which  are  discussed  below. 

The  first  way  to  produce  figure  6 is  to  select  a_.  , b . , and  c_.  , 
i = 1,  2,  3,  4,  which  satisfy  the  following  conditions. 

(1)  an  t a12,  a13  = a14,  a^  i a13,  a12  t al3 

(2)  a21  - a22  = a?3  = a24 


(3)  1’13  + h12),  b14  i b13 


(4)  b21  = b22 


b23  = b24 


(17) 


(5)  c13  i 2(cn  + c12) 


(h)  c21  - c22  = c23 


'24 


With  the  A^s  defined  in  this  way,  no  individual  two- flat  provides 

estimates  of  any  of  the  main  effects.  Thus  a design  of  this  type  could 

be  resolution  IV  only  if  the  treatment  combinations  in  the  union 
4 

T = U (t  | A.t  = c.)  work  together  to  separate  main  effects  from  each 
i = l 

other  and  from  interaction  effects.  Unfortunately,  the  rigidity  of 
conditions  (2),  (4),  and  (6)  in  (17)  ensures  that  some  main  effects 
and  interaction  effects  will  be  aliased  in  identical  ways  on  all  four 
two- flats.  Suppose,  for  example,  that  the  a.,  b_. , and  c^,  i = 1,  2,  3,  4, 
are  chosen  according  to  (17)  to  give  systems  of  equations  (18). 
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(1  ) 


(3) 


[l  0 0 1 

|o  1 1 1 


t = 


10  2 1 

t = 

Y 

uii 

l 

(2) 


(4) 


10  11 
0 111 


10  2 2 
0 111 


t = 


t = 


(18) 


It  follows,  for  instance,  that  on  every  flat  the  main  effects  of  Factor 
2,  say  F^,,  are  aliased  with  the  F^F,  interaction  of  Factors  3 and  4 
according  to  the  correspondence  of  levels  given  in  (l9). 


r2(t2) 

0 

1 

2 


F F ( t.  + t ) 
3 4 3 4 

1 

0 

2 


19) 


That  is,  for  every  t_  - (t^,  , t3,  ) ' in  T the  following  are 

satisfied : 

Whenever  t0  = 0,  t3  + t = 1. 

Whenever  t?  = 1,  tg  + t = 0. 

Whenever  = 2,  t3  + t^  = 2. 

Thus  the  design  determined  by  systems  of  equations  (18)  does  not 

separate  the  effects  of  F^  from  the  interaction  of  F3  and  F^. 

In  general,  since  the  preceding  example  illustrates  a consistent 
occurrence,  designs  specified  by  (17)  cannot  be  resolution  IV. 

The  second  construction  that  yields  Figure  6 can  be  achieved  in  a 
number  of  ways,  a typical  one  being  specified  by  the  conditions  (20). 


(l) 

all  * 

ai2’  al3 

a 1 4 ’ ail  * 

ai3’  ai2  * ai3 

(2) 

a21  * 

a22 ’ a23  = 

a24’  a21  ^ 

a23’  a22  ^ a23 

(3) 

b22 

2b,,  t b_, 
11  21 

+ b12’  b13 

= 2b11  t 2b1?, 

2b12’ 

b14  ^ bl3’ 

b24  = 2bll 

+ b21  + b14 

(20) 


iiiii  »ii 


22 

= 2cn 

+ C21  + 

'12’ 

C 1 3 

13 

’ C24  Z 

?Cll  + 

C21  + 

' 14 

2 3 


:11  + ?C12’  C23  ?cll  + C21  + 


Designs  of  the  type  (20)  r'er.emble  those  of  type  (17)  in  that  1 he 

individual  flat"  do  not  guarantee  est imability  of  any  of  the  main  effects 
Moreover,  in  any  specific  instance  it  is  easy  to  see  that  t he  flats  do 
not  work  together  to  create  a resolution  IV  design.  Tor  example,  if  the 
d-L’  b^  , and  c^.  , i - 1,  2,  3,  4,  are  chosen  according  to  (20)  to  give 
systems  of  equations 


(1) 


(3) 


’l 

0 

0 

ll 

[ll 

(2) 

*1 

0 

1 

2' 

2l 

0 

1 

0 

2j 

t = 

.2. 

1 

1 

0j 

1 = 

L J 

1 

1 

0 

2 

? 

t = 

I1' 

(4) 

1 

o 

2 

r 

1 = 

°1 

.0 

1 

2 

lj 

L2. 

.0 

1 

2 

2_ 

Ll  J 

(21) 


it  follows  that  main  effects 


are  aliased  with  interaction  effects 


by  the  same  correspondence  on  all  four  planes.  Thus  Design  (21), 
like  Design  (18)  cannot  be  resolution  IV. 

A second  design  constructed  according  to  (20)  consists  of 
solutions  to  systems  of  equations  (22). 


(3) 


[lOOl] 

[0112. 

i°n 

Lo. 

j (2) 

r 1 01  7 
Lo  1 2 0. 

t = ! 

0 

0. 

1 0 2 0] 
.0  1 0 lj 

It  = 

[1] 

(4) 

1021 

.0112. 

t = 

r 3 

0 

.0. 

(22) 


The  aliasing  patterns  determined  by  these  four  flats  are  not  as  easy  to 

analyze  as  those  of  (21).  However,  to  cite  one  example,  on  the  first 

2 

flat  the  main  effects  for  are  aliased  with  interaction  l'^f^  according 


to 


0 0 
1 -<->■  2 

2 ■<— > 1 


while  on  the  second  flat  the  main  effects  of  F,,  are  aliased  with  F ^ F^ 
by  tho  same  correspondence.  This  occurrence,  along  with  other  similar 
correspondences,  suggested  that  Design  (22)  is  not  resolution  IV,  a 
result  that  has  been  verified  by  direct  computer  check. 

in  summary,  although  the  discussion  of  the  preceding  paragraphs 

4 

does  not  constitute  a proof  that  no  21  run  3 resolution  IV  design  can 
be  constructed  from  four  intersecting  two-flats,  it  doer,  indicate  that 
this  is  probably  the  case.  Moreover,  it  surely  points  out  the  need  for 
a means  of  analyzing  how  the  choice  of  A.  and  c . , i = 1,  ....  n,  affect 


n 
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4 . 1_.  t_i in.ih  1 1 1 1 y oi  r t t ■ -•i  1 

The  example  s md  discussion  of  the  previous  section  illustrate 

n 

nearly  how  the  est  imabi 1 ity  oi  effects  from  de:  ign  T = U (t_|A.^=c_. } 

i=l 

:■  pei  : nol  only  On  the  individual  A.s,  but  also  or.  the  tructure  oi 
inter  - t ion  resulting  from  choice  of  c. , i : 1,  2,  ....  n.  If  the  A.s 
are  : pec  if i < • i by  Theorem  2 of  Section  2,  it  is  of  course, 

poo  i'  !-■  to  estimate  all  main  effects  using  just  one  flat  at  a time. 

H wev  r,  tl  t imal  ility  of  p that  occurs  if  c,  = c^  = ...  = = £ 

: ided  not  by  any  individual  flat,  but  rather  by  the  joint  appli- 

cation o’  all  treatment  combinations  in  a union  of  n flats.  Similarly, 
the  t , . vet  • • de:  igns  summarized  in  Example  8 of  Section  3 are 

all  resolution  IV  although  estimability  of  certain  main  effects  is  not 
guaranteed  by  considering  the  flats  one  at  a time. 

The  present  section  introduces  an  approach  to  the  analysis  of  t lie 
aliasing  structure  for  designs  formed  as  the  union  of  a number  of  flats. 
The  discussion  includes  examples  which  illustrate  use  of  the  basic 
procedures . 

For  the  3n  experiment  consider  again  the  general  design 

r 

T = U (t  | A.t  = c . },  (23) 

i=l  1 

where  each  A.  is  an  m.  x n matrix  of  rank  m..  for  each  i,  i r 1.  ..  *r, 
1 1 1 

the  model  equat  ions  that  correspond  to  the  N.  - ' t r e.itmeiit  coinhini- 

tion  specified  by  A.t_  = c_.  are 

VA  V.  ) = X.g 

— L l — 


(24  ) 


26 


where  X.  is  an  N . x v matrix  of  known  constants,  and  0 is  a v x 1 
vector  of  non-negligible  unknown  parameters.  Thus  the  model  for  the  N 
observation.;  corresponding  to  treatment  combinations  in  (21)  can  !><■ 
partitioned  as 

0 = XB.  (25) 


for  every  i,  i = 1,  r,  the  coefficient  matrix  A.  partition 

the  parameters  in  0 into  a number  of  alias  sets,  sav  5..,  ....  S. 

— 1 1 i u ; 

4 

For  example,  m a 3 experiment  involving  factors  F,  , F?,  F.(,  and  F^ 
the  coefficient  matrix 


0 1 

1] 

Lo 

1 2 ■ 

J 

produces 

the  five  alias  sets  = 

{p}. 

S2  = {Fr  F2F3’  F3F4’  F2F4}' 

“3  " ^ 2 ’ 

br2V  rir4-  r3r4b  S»  - 

ir3’ 

FlV  FoF4  ’ F!F2>’  S5 

= ^4 » F^F^,  F'prg },  where  p has  one  degree  of  freedom  and  everv 

other  effect  has  two  degrees  of  freedom  corresponding  to  the  linear  and 
quadratic  contrasts  that  were  discussed  in  Section  1. 

Suppose  now  that  E is  an  effect  of  interest  and  that  estimability 
of  the  linear  (E^  ) and  quadratic  ( E^ ) components  of  E is  not  guaranteed 
by  any  one  flat  in  (3.23).  Let  Si  = {E  = E.,,  E.„ E.„  ) be  the 

t i 1 ' i2  it; 

alias  set  including  E that  is  generated  by  coefficient  matrix  A., 

i = 1,  ...,  r,  and  assume  that  each  E^.  . , j = 1,  ....  f . , i = 1,  ....  r, 

carries  two  degrees  of  freedom.  Then  if  E..  and  E...  are  any  two 

ill]' 


■ 

Y_ 

r 

x. 

-1 

i 

Y„ 

x„ 

-2 

2 

Y 

X 

— r 

r 
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member:',  of  SV , it  is  possible  to  transform  the  levels  of  E.  . to  the 
i 1J 

levels  of  E ^ , by  applying  one  of  the  permutations  in  (26). 


E.  . 

1 I 

0 

1 

2 


fill 

0 

1 

2 


(01): 


E.  . 

n 

0 

1 

2 


(012): 


(02): 


( 26) 


(021) 


(12): 


Moreover,  if  X.’  represents  the  submatrix  of  X.  consisting  of  the 
i r i 

columns  that  correspond  to  linear  and  quadratic  components  for 


E . 


U 


= 1,  . £.,  it  is  clear  that  X.  is  of  rank  two  since  exactly 


two  degrees  of  freedom  are  associated  with  each  effect  in  S.. 

Suppose  specifically  that  linear  and  quadratic  contrasts  are  specified 
by  the  orthogonal  polynomial  coefficients 


Then  if  E.j  and  E..,  in  si'  are  related  by  permutation  g of  (26),  and 
if  Xf:  ...  consists  of  the  four  columns  of  X;  associated  with  E.  . and 

i;n  1 n 

E..,,  one  finds  that 
ij 


X.;jj,  Pg  - 0, 


where  Pg  is  determined  by  permutation  g in  accordance  with  (27). 


& 


e 


(012) 


(021) 


Pg 

2 0 
0 2 
-2  0 
0 -2 

2 0 
0 2 
1 2 
-1  1 

2 0 
0 2 
1 -3 
1 1 
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(01 ) 


(02) 


(12) 


Pg 

2 0 
0 2 
-1  -3 
-1  1 

2 n 
0 2 
2 0 
0 -2 

2 n 
0 2 
-1  3 
1 1 


( 27) 


It  is  instructive  to  consider  first  the  special  case  for  which 
A.  = B,  a constant  m x n coefficient  matrix  independent  of  i.  Then 
design  (23)  becomes 

r 

T = U (t_  | Bt_  = d_.  } , (28) 

i = l 1 

each  flat  consists  of  N.  = 3n-m  treatment  combinations,  and  the 

resulting  union  of  parallel  (n  - m)-flats  includes  a total  N = (rX3n~m) 

r r 

treatment  combinations.  Moreover,  S.  = S = {E  = EJ?  E?,  ...,  } 

is  t he  same  for  all  r flats. 

Consider  for  each  i,  i = 1,  ....  r,  the  null  space  of  X*.  . That 
is,  given  that  includes  £ effects,  each  with  two  degrees  of  freedom, 
consider  the  space  of  all  x 1 vectors  w such  that  X.w  = 0. 

I - 

Theorems  3 and  A,  which  follow,  establish  a means  for  investigating. 

the  e.stimability  of  F and  E through  consideration  of  the  null  spacer; 

Li  Q 

for  X? , i = 1 , . . . , r. 


29 


£ 

Theorem  3 . Let  WV  be  a 2l  x 2(£  - 1)  matrix  of  full  column  rank  such 

E E 

that  XV  W.  = 0,  where  0 represents  an  N . x 2(1  - 1)  matrix  of  all  zeros. 

E E 

Then  the  column  space  of  W.,  denoted  C(W.‘),  is  the  same  as  the  nuil 

l l 

P VE 

space  of  X.. 

E . E 

Proof.  Eince  XV  is  of  rank  two,  if  w is  a solution  to  XVw  = 0_,  it  is 

possible  to  solve  for  two  components  of  w_  in  terms  of  2f  - 2 arbitrary 

p p 

components.  Thus  the  null  space  of  XV  has  dimension  2 £ - 2.  If  WV 
is  of  full  column  rank  and  such  that  X^V  W‘V  = 0,  the  columns  of  WE 

ii  i 

p 

clearly  form  a basis  for  the  null  space  of  XV.  The  proof  is  complete. 


Theorem  4 . The  linear  and  quadratic  components  for  E are  not  esti- 
mable from  the  runs  of  design  (28)  it  and  only  if  there  exists  a 

r p 

vector  w e (*)  C(W.)  such  that  the  comDon‘*nts  of  w which  correspond  to 

i = l 1 

E.  and  En  are  nonzero. 

L v 

Proof . In  terms  of  model  (25)  any  component  of  8 is  estimable  if 
and  only  if  the  corresponding  column  of  X is  linearly  independent  of 


all  remaining  columns  of  X.  Since  the  structure  of  design  (28) 
guarantees  that  E.  and  E can  be  aliased  only  with  effects  in  S1  , it 

L (J 

suffices  to  consider  only  the  submatrix  of  X consisting  of  columns  for 
the  linear  and  quadratic  components  of  effects  in  S'.  Partition  this 


X 


E 


submatrix  as 


30 


r 

and  lot  the  first  two  columns  of  X“  correspond  respectively  to  E and 

E . How  E is  not  estimable  if  and  only  if  there  exists  a 2f  x 1 
x L 

vector  w = [w^ , w?,  w^]'  such  that  x'w  = £ and  in  which  w i 0. 

E 

Clearly  such  a w must  he  in  00(W.').  A similar  argument  holds  for  E(, , 
and  the  proof  is  complete. 

E 

In  practice  the  matrix  W.  which  provides  the  basis  for  the  null 

space  of  xl  can  be  formed  in  any  of  a number  of  ways.  Two  of  the  most 

useful  form>  are  specified  bv  Definitions  2 and  3. 

E E 

Definition  2.  W.,  a basis  for  the  null  space  of  X.,  is  said  to  be  in 

l i 

standard  form  with  respect  to  effect  E,  , where  E,  may  or  may  not  be 

k k 

equal  to  E,  if  wl‘  is  written  as 


"2  0 
0 2 

2 0 
0 2 

2 0 
0 2 

2 0 
0 2 

. . . 

2 0 
0 2 

Zkl 

0 

0 

Zk2 

0 

0 

0 1 

• i 

0 

• 

• 

• 1 

. 

• 

° r ° 

csr* 

0 

* 

0 

• 

• 

« 

.0 

0 

0 

0 

• 

0 

zk£  J 

where  each  0 block  is  a 2 x 2 matrix  of  zeros  and  each  2,  . is  deter- 

kl 

mined  from  (27)  according  to  the  permutation  required  to  change  the 

levels  of  E,  to  the  levels  of  E.,  j = 1,  l,  j ^k. 

k 1 

E E 

Da!  in  i r iori  3 . W.,  a bas  is  for  the  null  space  of  X.,  is  said  to  be  in 

p 

diagonal  form  if  Wj‘  is  written  as 


31 


where  for  some  ordering  of  the  effects  in  S ,7,,  , „ is  determined  from 

k,  k+1 

( 27  ) according  to  the  permutation  required  to  change  the  levels  of 

C,  to  the  levels  of  C,  ,,k=l,  ...,£-  1. 
k k+1 

3 

Example  9.  For  the  3 experiment  involving  factors  F . F„ , and  F',, 

7 1 * 

consider  design  T = U (_t_  | Bj_  = cl_.  },  where 


and  dj  = (0,  0)’,  d?  = (0,  1)’,  d_3  = (0,  2)’,  d^  = (1,  0)',  d,.  = (2,0)’, 
d^  = (1,  1)’,  and  d_7  = (2,  2)'.  The  resulting  design  T consists  of  the 
21  treatment  combinations  shown  in  (29),  where  spacing  indicates  the 
seven  individual  flats. 


‘012 

012 

012 

01? 

012 

012 

01?j 

T = 

01  2 

012 

01? 

120 

201 

120 

201 

(29) 

012 

120 

201 

01? 

012 

120 

201 J 

Coefficient  matrix  1 

$ producer. 

alia 

n set 

o/l 

= sr? 

. 

— o 

3 r r 

1 iq,  F2,  r3,  F,r2,  r,r3,  r,r3>  .ho  w'.i  * wb  , ^3  = Wj> 

i r 1,  7,  eacli  written  in  diagonal  form,  are  as  given  in  (30). 

Unspecified  entries  in  each  W.  are  all  equal  to  0. 


i. 


2 0 

0 2 

2 0 2 0 

0-202 
-2  0 


0-2  0 


- 1 3 2 0 

110  2 


-2  0 2 0 
0-202 
-2  0 
0 -2 


110 

2 


0-2  0 


-110  2 

2 0 


-2  0 
0 -2 


2 0 2 0 

0-2  0 ->2 

13  2 0 


2 0 
0 2 

-2 

0 

2 

0 

0 

-2 

0 

2 

1 

3 2 

0 

-1 

1 0 

2 

-1 

2 

0 

-1 

1 

0 

1 

2 

j q n 

i 

-1 

1 0 2 
-2  0 

A 0 

~ 

W„ 

u - / 

2 

~ 2 

0 

— 

0 

2 

1 

3 

2 

0 

-1 

1 

0 

2 

1 

-3  2 

0 

1 

1 0 

2 

-1 

-3 

2 

0 

1 

1 

0 

2 

1 

-3  2 0 

1 

1 0 2 
1 3 

i i 

- 

w 

— l l 

4 

2 

0 

- 

0 

2 

1 

3 

2 

0 

-1 

1 

0 

2 

-2 

0 2 

0 

0 

-2  0 

2 

2 

0 

o 

c. 

0 

0 

-2 

0 

/ 

33 


2 0 
0 2 

1-320 
110  2 
- 2 0 2 0 
0-202 

-1-3  2 0 

-110  2 


0 2 0 


0-202 
1 -3 


By  Theorem  4 for  k = 1,  2,  or  3 the  linear  and  quadratic  effects  of 
factor  F^  are  not  estimable  if  and  only  if  there  exists  a vector 
w = [w^ , w_,,  w j 2 J ' with  nonzero  components  corresponding  to  ^ 

and  F,  _ such  that  w can  be  written  as  some  linear  combination  of  the 

k i" 

columns  of  W^,  say  W^a.,  for  each  i,  i = 1,  7.  Thus  to  prove  the 

estimability  of  F r and  F.  _ it  suffices  to  show  that 
K 5 Lj  k j Q 


w ~ = W2-2  = • • • = w73.7 


'll’  fi12' 


•’  “l,10]'*  - = Vl 


forces  the  components  of  w which  correspond  to  F 
zero.  Now  if  = [a^, 
implies  that 

rr« 

?ai2 

■3ail  “ 3al2  + Jai3 

”ll  - 3nl  2 + 

-3all  ‘ 3a12  + fl5 

ot^  - -’nt.o  + 2a,, . 

3all 

3aH 
3all 

all 


, . and  T.  . to  be 
k,L  k ,0 


wrA>  = w7“7 


+ 

s;  - 

- iu 
3(1 15 

+ 

3v, 

+ 

2a 

“ 

3al  2 + 

“l5 

+ 

« 6 

2a 

+ 

3rll  2 ' 

3al  5 

+ 

3«1B 

2a 

- 

3al  2 + 

a15 

+ 

3a16 

~ 

2a 

17 

18 
19 

1,10 


an>  1 


i 


34 


—7 


2 a 
2a 
-3a 


<l> 


-3a 


-3a 


-3a 


11 

12 

11 

11 

11 

11 

11 

‘11 

11 


L u 


+ 3a 

- 3a 
+ 3a 

- 3a 

- 3 a 
+ 3a 

- 3a 
+ 3a 


2a 
2a 
2 a 
2a 
fia 
2a 
6a 
2a 


13 

14 

15 
15 
lb 

15 

16 
15 


2a 

2u 

2a 

2a 


17 

18 
19 


1,10 


Thus  wu  = -2«ig  = 3au 


3al2  + 6a!6  + a!9  + 3al,10  = -6a17  - 3a^ 


15 


+ al9  - 3al,10  3nd  W12  = -2al,10 


11 


3ct12  + ?a15  - V>  + al  ,10 


= -2a,,  + a, c + 3a, _ + a,_  + 2a,  so  the  equations  of  (31)  must  he 
11  15  16  19  1,10 


satisfied . 
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Theorem  5 . Let  the  3n  design  T = f]  {t  | Bt  = d.},  where  B is 


(n  - 1)  x n matrix  of  the  form 


and  d.  = (d^,  d^,  d.  n_1)',  i = 1,  r.  for  any  k, 

k = 2,  n,  the  linear  and  quadratic  effects  of  the  kth  factor  are 

estimable  from  the  runs  of  T if  the  set  of  all  d_.  , i = 1,  r, 

includes  a subset  of  size  three,  say  d , d and  d , in  which  d = d 

J — s — u — v sy  uy 

= d for  y = 1,  ...,k-l,k+l,  . n - 1,  while  the  kth  component 

ranges  over  the  three  values  in  GF(3). 

Proof . Consider  the  set  of  treatment  combinations 

= U { t_  | Bt_  = d^.  } , where  B is  defined  by  (23)  . Since  d , 
.i=s,u,v 

, and  dv  are  equal  except  for  the  kth  component,  which  ranges  over 
all  three  possible  values,  the  levels  of  Factor  k are  not  associated 
with  the  levels  of  any  other  effect  by  one  of  the  permutations  in  (261. 
Thus  given  the  model  E(Y^)  = where  is  the  vector  of  random 

observations  correspond ing  to  treatment  combinations  in  T^4 it  is  not 
possible  to  involve  the  columns  associated  with  linear  and  quadratic 
F^  in  a linear  combination  equal  to  £.  Thus  in  the  model  for  the  full 
design  T the  columns  of  X associated  with  F are  linearly  independent 

K 

of  all  remaining  columns  of  X,  so  the  linear  and  quadratic  effects  of 


F,  are  estimable.  The  proof  is  complete. 


36 


Suppose  that  the  conditions  of  Theorem  5 are  satisfied  for  each 
k,  k = 2,  . . . , n.  If  the  condition  is  satisfied  for  k = 1 after  taking 
the  appropriate  linear  combination  of  the  rows  of  B and  the  elements  of 
d_.,  i=l,...,  r then  the  design  is  resolution  IV.  Example  10  illus- 

4 

trates  how  this  approach  can  be  used  to  consider  estimability  in  a 3 
foldover  design  of  the  type  discussed  previously  in  Section  3. 

4 

Example  10.  A 3 foldover  design  for  a factorial  design  involving 

factors  , F F^,  and  F^  was  given  in  Section  3 as  a union  of  four 

intersecting  flats.  This  design  can  be  expressed  in  the  format  of 

9 

Theorem  5 as  T = U { t_  | Bt_  = d_. } , where 
i=l 


2 10  0 

B = 2 0 1 0 

2 0 0 1 


and  the  d.,  i = 1,  ...»  9,  are  the  columns  of  the  matrix  (33). 


000001212 

000120012 

012000012 


[— 1 ’ -2’  - ’ -9-1 


Estimability  of  the  main  effects  of  F? , Fg,  and  F4  is  guaranteed  by 
Theorem  5 as  follows:  F ^ by  d^ , dg  and  d_7;  F^  by  d_^  , d^  , and  d?; 

F3  by  d_^ , d_2,  and  d3>  Since  elementary  row  transformations  will  trans- 


r 

2 

1 

0 

0 

0 

1 

*1 

0 

= 

2 

0 

1 

0 

0 

1 

2 

2 

0 

0 

1 

0 

1 

2 

"l 

2 

0 

0 

2 

1 

0 

0 

2 

1 

0 

0 

0 

0 

_0 

2 

0 

1 

0 

0 

0_ 

the  main  effects  of  F^  are  estimable  also.  Thus  the  design  is  resolu- 
tion IV. 


Theorem  6.  Suppose  that  d ^ of  (28)  are  selected  so  that 
for  some  in  S*'  = (Ej  , E^ , ...,  E^>}  the  following  conditions  are 
sat isf  ied. 

1.  For  each  j i k , the  permutation  of  (26)  that  relates  E,  to  E.  is 

K 1 

the  same  for  each  of  d,  , d„,  . . . , d . 

—1—2  — r 

2.  The  list  of  r permutations  relating  E^  to  E for  d ^ , ...,  d^, 
respectively,  includes  at  least  two  distinct  permutations. 

Then  design  (28)  ensures  the  estiroabili ty  of  the  linear  and  quadratic 
effects  of  E. 

E 

Proof.  If  necessary,  reorder  the  elements  of  S so  that  E^  = E^  and 

E = E^-  For  each  i,  .i  = 1,  ...,  r,  write  a basis  VE  for  the  null 
E 

space  of  in  standard  form  with  respect  to  E^  = E^ . That  is. 


where  V..,  j = 2,  3,  ...,  ( is  determined  by 
permutation  relationship  between  E^  - E^  and 


(27)  according  to  the 
r 

E. . Mow  w e H C(W. ) 

? - 1 


I 


Since  V.a.v  = [>5,  w^]'  for  all  i,  a.y  = V^[ 

all  i,  i = 1,  r.  Thus 


effects  of  L are  estimable.  The  proof  is  complete. 


Fxam:  !»■  li. 
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consider  T 


For  the  3 experiment  involving  factors  , 
3 

= U { t | Bt  = d . } , where 

i=l  - - -1 


f2,  r3,  and 


B 


[l  0 1 \ 

[0121 


and  d_^  = (0,  0)',  d,,  = (1,  0)',  d_3  = (2,  0)'.  Coefficient  matrix  B 

p cp 

generates  the  alias  et  S 1 = { F^  , F^F^,  F^F^,  F^F^},  arid 

021112200  102220011  210001122 

012012012  012012012  012012012  . 

012100221  012100221  012100221  '• 

000121212  000121212  000121212 


where  tne  spacing  indicates  the  treatment  combinations  on  the  indivi- 
dual flats.  From  (34),  one  sees  that  the  permutation  relating  F^F^ 

to  Fj  is  e on  flat  1,  (012)  on  Flat  2,  and  (021)  on  flat  3.  Since 

2 

F^Fj  is  related  to  F F|(  and  to  F^F^  by  the  (12)  permutation  for  all 
runs  in  T,  Theorem  0 guarantees  that  T provides  estimates  for  the 
linear  and  quadratic  effects  of  F^ . 

The  preceding  discussion  of  the  aliasing  structure  for  designs 
formed  as  the  union  of  r parallel  flats  extends  readily  to  designs 
specified  by  (23)  as  a union  of  intersecting  or  non  intersect ing 
flats.  In  this  more  general  case,  however,  the  alias  set  sl‘ 


= {Eil  = K’  E'i2’ 


X!:  is  taken  to  be 


* E . / 

i^i  that  includes  effect  F.  depends  on  A.,  and 
the  submatrix  of  X.  consisting  of  all  columns  that 


l i 

correspond  to  effects  in 


,Q  . 


r 

U 

i=l 


fan 

J = 1 1 


l 


Jk 


4U 


where  0 - U S*?.  That  is,  consists  of  all  effects  that  are  aliased 
i = l 1 

with  an  effect  in  Q on  any  flat,  i = 1,  ...,  r.  It  is,  of  course, 
permissible  to  eliminate  from  Sf  any  effect:  that  is  known  to  !>•• 
estimable  from  the  treatment  combinations  on  one  of  the  individual 


flats.  For  each  i,  i = 1,  r,  the  W.  matrix  that  supplies  a basis 

F, 

for  the  null  space  of  is  con. posed  of  one  or  more  matrices  of  the 

type  specified  in  Definitions  2 and  3,  depending  on  how  many  of 

Q 

the  alias  sets  generated  by  A.  are  included  in  S . The  notions  of 

this  paragraph  are  illustrated  in  Example  12. . 

Example  12.  Design  3 of  Example  8,  Section  3,  is 
4 

T = U (t  I A.t  = ol,  where 
i = l “ ~ 


4 • IS 


2 

1 

0 

A - 

'2 

0 

1 

0 

2 

0 

lj  ’ 

2 

_7 

0 

0 

1. 

1 

0 

0 

A = 

'l 

0 

1 

1" 

0 

0 

1 J ’ 

4 

0 

1 

1 

1, 

all  over  GF(3).  Theorem  1 guarantees  the  estimability  of  the  main 
effects  of  Factors  F^,  F , and  F^  individually  from  flats  1,  2,  and  3, 
respectively.  Estimability  of  the  linear  and  quadratic  components  of 
F[( , on  the  other  hand,  does  not  follow  immediately  from  consideration 


of  any  one  flat.  In 


fact,  s[4  = (F4,  F2,  F3,  F2Fv  F2F4,  F3r4}, 


S24  = (r4,  r1 , F3,  Fj f3 , riF4,  r3r4>,  s3“  = lr4,r1,r,.rJr,,,r1r4,r.,r4l, 

and  s|(4  = (F4,  FjFg,  F^F 3 } . Thus  0 = U S1.’4  = {F  , F„,  F Fl( , F V , 

5 

F^Fg,  FjF4,  FgFg , F?F4,  FgF4 } . Since  F^ , F?,  and  F^  are  known  to  he 
estimable  from  flats  1,  2,  and  3 individually,  take 

^4’  fl^3'  I'lt4’  ^3’  1 4’  1 3f4^‘ 


1 


41 


r1 

The  alias  sets  needed  to  construct  W.  , i = 1,  2,  3,  4 and  thereby 


complete  the  analysis  are 


f’lr2 


= {F1F2’  F1F3>  FiV»  S2 


rii’2 


{ 1 ! 1 2 ' F2F3’  *2*4^' 


C-F21'3  - 


21 3’  1 11 3’  1 31 4 f ’ 


(-rlF2_ 


= {F1F2’F3V’ 


and  SFl‘4  = {I’ir4,  The  w!4 , i = 1,  2,  3,  4,  which  are  displayed 

in  (35),  are  composites  of  matrices  constructed  in  the  standard  form 
specified  by  Definition  2.  The  unspecified  entries  in  each  wl4  are 
all  equal  to  0. 
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the  equations  of  (37)  must  be  satisfied. 
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Elementary  row  operations  will  reduce  (37  ) 


to 


10  0 
0 0 1-1 
0-1-2  1 
0 11 
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or 
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15 
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which  implies  that  a^  =0.  A similar  consideration  of  wu , w^,  and 


implies  that  a^ 0 = 0.  Thus  w1  = -a1 1 + 3a 


12 


1 


11 


12 


0 and  w„ 


ail+f*12 


= 0,  so  the  linear  and  quadratic  effects  of  F are  estimable. 
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